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Abstract. The rich-get-richer mechanism (agents increase their "wealth" randomly 
at a rate proportional to their holdings) is often invoked to explain the Pareto power- 
law distribution observed in many physical situations, such as the degree distribution 
of growing scale free nets. We use two different analytical approaches, as well as 
numerical simulations, to study the case where the number of agents is fixed and finite 
(but large), and the rich-get-richer mechanism is invoked a fraction r of the time (the 
remainder of the time wealth is disbursed by a homogeneous process). At short times, 
we recover the Pareto law observed for an unbounded number of agents. In later times, 
the (moving) distribution can be scaled to reveal a phase transition with a Gaussian 
asymptotic form for r < i, and a Pareto- like tail (on the positive side) and a novel 
stretched exponential decay (on the negative side) for r > \. 
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1. Introduction 

The rich-get-richer mechanism, or preferential attachment, is often invoked to explain 
the ubiquity of Pareto distributions in complex networks and in natural and man-made 
phenomena [HE], [3]. For example, in the Barabasi- Albert model for complex networks [4] 
a power-law tail in the degree distribution of the nodes, P(k) ~ k~ 3 , arises when new 
nodes are attached to the net, one at a time, to one of the existing nodes with a 
probability proportional to their degree. Crucial to the emergence of this power-law 
tail, however, is also the fact that the net grows with time (new nodes are attached) - 
preferential attachment alone is not the only necessary ingredient. 

In this paper we study the consequences of preferential attachment without growth, 
that is, when "wealth" is disbursed to a finite group of A agents by the rich-get-richer 
mechanism. More precisely, wealth is incremented by the rich-get-richer mechanism only 
a fraction r of the time; the remainder of the time (a fraction 1 — r) wealth is disbursed 
homogeneously at random. For an infinite number of agents, 1 — r is simply the rate at 
which new agents are introduced (are given their first unit of wealth), and there results a 
Pareto distribution of wealth, P(k) ~ k~ x , with A = 1 + 1/r. When A is finite, however, 
the homogeneous process may allot wealth to agents that had already been introduced. 
The case when r = and wealth is disbursed only by the homogeneous process results 
in a Gaussian distribution of wealth — quite different from the power-law distribution 
found in the limit of A — > oo. The two extremes, of A — > oo, and A < oo and r = 0, 
dominate the general case, pulling the system toward opposite poles and resulting in 
a kinetic phase transition: the long-time asymptotic distribution of wealth is Gaussian 
for r < |, but has a power-law tail for r > |. In all cases the wealth distribution can be 
written in scaling form, and the shape of the scaling function is characterized by various 
r-dependent exponents that we obtain analytically. 

The rest of this paper is organized as follows. In section [2] we present the model 
and analyze the two special limits that color its behavior. This section serves also to 
illustrate the two mathematical techniques that we employ for the model's analysis: the 
master equation (for the number of agents with a given wealth) and the rate equation 
approach (called also mean-field by other researchers). In the latter approach one uses 
a rate equation to compute the expected wealth of a particular agent as a function of 
the time of its introduction. The wealth distribution may then be obtained from the 
distribution of the introduction times. The general case is treated in section [3j We 
derive several analytical results and compare to computer simulations, as well as discuss 
the range of validity of the rate equation approach. We conclude with a summary of 
our findings and a discussion in section HI 

2. The model 

We analyze the following toy model of wealth allotment. A single wealth unit is given per 
unit time to a random member out of A agents. With probability r the agent is chosen 



Phase transition in rich- get-richer mechanism 



3 



by the rich-get-richer mechanism, proportional to its current wealth; with probability 
1 — r the agent is chosen uniformly at random. Thus r is a tunable parameter controlling 
the relative strengths of these two processes. We assume that initially each agent has 
zero wealth. Let us start by exploring the two limits of A — > oo and r — > 0. These 
exactly solvable extremes yield markedly different behavior and underlie the essential 
physics of the system when A and r are finite: for r < | the uniform random process 
dominates, while for r > | the rich-get-richer process gets the upper hand. 



2.1. The limit of A — > oo 

For A — > oo, agents with positive wealth form a set of measure zero for any finite time. 
Thus choosing an agent uniformly at random is equivalent to the "birth" of a new agent i 
with wealth fcj = 1. We analyze this case with a master equation approach, following 
the techniques and notation of [5]. The number of agents with wealth k > at time t, 
N k (t), obeys the master equation: 

±-N k = (1 - r)6 kA + „ r — \(k - l)JV fc _! - kN k ] . (1) 

Since one wealth unit is disbursed per unit time, Yl k kN k = t. This limit is characterized 
by both growth and preferential attachment, hence we expect a power-law distribution 
of wealth. This master equation has been studied before [S] and indeed it is easy to 
show that the actual distribution of wealth is 

P{k) =N k (t)/N(t) ~ k- x - A = l + i. (2) 

One can also derive this distribution using a rate equation for the average wealth 
of agent i: 

jk.it) = r -h(t), (3) 

since J? . kj = t. Because A — ► oo, the rate at which agent i is selected by the 
homogeneous random process is zero, while the rate of selection by the rich-get-richer 
mechanism is kij kj = k./t. Following the analysis in [6], one can show that 

P(k) = -k~ l ~ l/T , k>l, (4) 
r 

i.e., a Pareto distribution with the same power-law tail. Note that this distribution is 
properly normalized (taking k to be a continuous variable) and that its first moment 
agrees with that from the master equation. 

While the results for both the master equation and the rate equation approach are 
equivalent for A — > oo, we will show that this is no longer true for A < oo. 

2.2. The limit of r -> 

When the number of agents A is finite, the N k obey the normalization condition 

oo 

A, (5) 

k=0 
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where now we include in the counting agents with zero wealth (k = 0), and the 
distribution of wealth is P(k,t) = Nf.(t)/A. The mean wealth per agent increases 
linearly with time: 

<*> = J ■ (6) 

Consider the limit of r — > 0, where wealth is disbursed only by the homogeneous random 
process. The corresponding master equation is 

with initial and boundary conditions Nk(0) = A5k,o and N-i(t) = 0. This is a simple 
Poisson process, as confirmed by the solution of ([7]): 

N k (t) = A^^e-^ . (8) 
For i>i we apply the Sterling approximation to obtain the distribution 

P(M) , 1 ~rf '-A^m (9) 

Thus, P(k,t) has a power-law tail in the limit A — > oo (section 12.1ft . but is Gaussian 
when A is finite and r — > 0. 



3. Finite A and r 

Master equation approach 

In the general case of A < oo and r > the master equation for the process is 

j f N k = ^ (jV fc _! - JV fc ) + ^ ((fc - 1)^ - fciV fc ) . (10) 

The system is then simultaneously pulled toward the two different limiting behaviors 
analyzed in section [21 We will show that for r > | the rich-get-richer mechanism 
dominates the process and the wealth distribution develops a power-law tail (as for the 
limit of A — > oo), while for r < | the homogeneous random process dominates and the 
wealth distribution tends to a Gaussian (as for r — > 0). Because A is finite, (k) = t/A 
increases linearly with time. The width of the distribution of k around the average 
grows like t a , where the scaling exponent a = r for r > | and a — ^ for r < ^. At the 
transition point, r = |, the width scales as \Aln£. 

To see this, begin by approximating the discrete distribution N k (t) by its continuous 
counterpart, P(k,t). Expanding to first-order, equation f TTUT) now reads 
d 1 — r 9 r d 

and the method of characteristics yields the scaling solution 

= r«/ (^^) ; a = r. (12) 
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This, however, cannot be true for all values of r, as it disagrees with the distribution 
found for r = 0, where the scaling exponent is a = | instead of a = r = 0. The reason 
for this discrepancy is that in this case the Kramers- Moyal expansion |7J of ( flOl) must be 
carried out beyond the first order. Indeed, on substituting the scaling form (fT2|) in the 
master equation (with unspecified a), and carrying out the expansion to second-order, 
we find 

(a - r)t 2a f{x) + {a- r)t 2a xf{x) + ^/"(x) = , (13) 

where prime denotes differentiation with respect to x = (k — t/A)/t a , and we have 
omitted terms proportional to t a (these are negligible compared to t 2a , as t — > oo). If 
a > |, the term proportional to t can be neglected in the long-time limit, and ([TBI) is 
satisfied provided that a = r. Thus, the scaling form (fl2|) is valid only for r > ~. For 
r < ~, however, the second-order term in (1131) may not be ignored. The only non-trivial 
way to cancel out the time dependence is then to have t 2a = t. Thus, for r < ^ the 
scaling exponent is a — |. At the transition point, r = |, there is no way to get rid of 
the time dependence in ( {TBI with the scaling form ( fl2i) . Taking a cue from other phase 
transitions we guess a scaling form with logarithmic dependence: 

P (it )f ) = J_/f*^/!!V r = I. (14) 

v ' ; (tint)" 17 V (tlnt) a / 2 v ; 

On expanding the master equation with this scaling form the leading behavior in time 
cancels out, provided that the scaling exponent is a = h. The next largest terms (smaller 
by a In t factor), yield the equation 

f(x)+zf'(x) + jf'{z) = Q- t r=i, (15) 

where now x = (k — t/A)/\ / t Int. In all three cases (for r), expanding to third- or 
higher-order yields additional subdominant terms. From the largest subdominant term 
one can deduce how fast the system reaches the scaling regime: the transient dies off as 
£-(2r- 1) f or r > i ag £-1/2 f or r < |, and as (hit)" 1 for r = |. Thus at the transition 
point, r = |, there occurs a critical slowing down as the system creeps into the eventual 
scaling regime logarithmically slow. 

For r < | we can use ( TTBl to find out f(x) and show that the limiting form of the 
wealth distribution is Gaussian: 

as t — ► oo. The divergence of the width of this distribution as r — > | is reconciled with 
the fact that at the limit r = ~ the scaling parameter picks up a (diverging) logarithmic 
component. The scaling function is still Gaussian, as can be deduced from (fT5l) : 



1 ,, Jt-i/A) 21 



r<2, (16) 



2 tint 



t/Af 



(17) 
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For r > | equation (fT3l) yields a tautology and one is unable to determine f(x). It is 
possible, nevertheless, to infer the limiting behavior: 

.-l-l/r 



/(*) 



cxp 



.T — > OO, -j 

r > - . (1* 

-(l-r)(A r \x\)"-~ 



— oo, 2 



The limit for x — > oo follows from comparing the distribution P(k,t) for the case of 
A — > oo with /(x)|a^oo = f(kt~ r ). For x — > — oo, we observe that the density of agents 
with zero wealth decays as N ~ exp[— (1 — r)t/A], see equation (|2"T|) . and we compare 
to f(x)\k=o = /(— ^ 1 ~ r M), leading to the second line of ffT8l) . An alternative derivation 
is presented next, using the rate equation approach. 

3.2. Rate equation approach 

The rate equation for the wealth of agent i, in the general case, is 

^**(*) = — + , ( 19 ) 

with initial condition fcj(tj) = 1. The solution, 

wK 1 -!)©'-^ < 2 °» 

is monotonously decreasing in £j. 

The probability ^(i) that an agent has still zero wealth at time t satisfies the 
equation 

|*( ( ) = (21) 

so ^i(i) = exp[— (1 — r)t/A]. It follows that the probability that agent i has been 
introduced (gets its first unit of wealth) by time T, given that it has been introduced 
by time t, is 

X(T) = \~ e Z? t • (22) 
1 — e a 1 

Note that this has the limit T/t, as A —>■ oo, that is normally used for this case [6]. 

Finally, P{k,t) = —dx(T)/dk, where T(k) is the solution to h(T) = k. Since pi 
cannot be inverted analytically (other than for special values of r), we express P in 
parametric form: P(k(T),t) = -dx{T)/dk = -(dh/ dU\ u=T )- l dx{T) / dT , and k(T) is 
obtained by putting tj = T in fl20l) . The wealth distribution in parametric form is then 

a.(T)= (\-t) T -\ f(T) = -^—^-e-^ T , (23) 
y ' \ A) ' JV ; Ar + (l-r)T ' v ; 

where we have used the scaled expressions x = (k — t/A)/t r and / = t r P, taking the 
limit of t —>■ oo at the end (the fact that the limit exists and is finite confirms this 
scaling) . 

It is now easy to verify the asymptotic behavior (TT8I) . The limit x — > oo corresponds 
to T — > 0. In this limit, the second equation of ( 1231) gives / ~ x 1+r . But since 
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{k-t/A)/{t\ntfl 2 

Figure 1: (color online) Scaling of the wealth distribution in each of the two phases 
at r = | (top left), r = | (top right), and at the transition point r = ^ (bottom). 
Convergence to the scaling forms is rapid for r = | and r = | but logarithmically slow 
for r = | — note that in the latter case the data (for exponentially increasing times) 
is slowly creeping toward the Gaussian limit of ( JT7|) (solid line). The theoretical limit 
of (fT6l) (solid line) fits the case of r = \ perfectly, but the prediction (1231) from the rate 
equation approach (solid line) fits the case of r = | only qualitatively (besides agreeing 
with the overall scaling). 

T ~ x" 1 ^, from the first equation, we conclude that / ~ x _1_1 / r . The limit 
x — » — oo corresponds to T — >• oo. In this limit, the second equation of (1231) gives 
/ ~ exp[— (1 — r)T/A], while from the first equation x ~ — (l/A)T 1 ~ r . We conclude 
that / ~ exp[-(l - r^A^xl) 1 ^ 1 -^}. 

Clearly, the foregoing rate equation method does not apply to < r < |, for it 
fails to reproduce the appropriate scaling forms in this range. Thus the rate equation 
approach is viable only when the second-order in the Kramers-Moyal expansion of 
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the corresponding master equation may be neglected. In figure [TJ we show numerical 
simulations for r below, above, and at the transition point. The results confirm the 
scaling forms found analytically above. For r < \ convergence to the Gaussian pdf is 
relatively fast, while the critical slowing down at the transition point, r = ~, prevents us 
from attaining the analytical limit (TT71) in practice. For r > \ convergence to the scaling 
pdf is again quick, however the explicit form predicted by the rate equation approach is 
correct only qualitatively: we ascribe this to the fact that the second-order is implicitly 
missing in this approach. 

4. Summary and discussion 

In summary, we have studied a model of a finite number A of agents that accrue 
"wealth" by the rich-get-richer mechanism a fraction r of the time (wealth is disbursed 
homogeneously randomly the remainder 1 — r of the time). In the early time regime, 
or, equivalently, when A — > oo there results a Pareto distribution of the wealth k: 
P{k) ~ kr x , with A = 1 + 1/r. In the long time asymptotic limit, the system is attracted 
to one of two opposite poles, and there is a kinetic phase transition as a function of the 
parameter r. If r < |, the distribution tends to a Gaussian of width t/[(l — 2r)A]. If 
r > ~, the distribution keeps its power-law tail ~ k^ 1 ^ 1 ^ for large k. 

In all cases the wealth distribution tends to an asymptotic scaling form as a function 
of x = (k — (k))/w(t), where (k) = t/A is the average wealth amassed by an agent by 
the time t, and w(t) = t a is a measure of the width of the distribution. The exponent a 
undergoes a phase transition: a = | for r < |, and a = r for r > ^. At the transition 
point, r = |, there appear logarithmic corrections: w(t) = (tint) 1 / 2 . 

The scaling form of the wealth distribution f(x) = t r P in the regime r > ^ is 
characterized by two more exponents (in addition to the width exponent a = r): 
f(x) ~ x~ l ~ l l r for x — ► oo, and f(x) decays as a stretched-exponential, with power 
1/(1 — r), as x — > — oo. Finally, the approach to the eventual scaling form ~ t~ z is 
characterized by a fourth exponent: z = \ for r < |, and z = 2r — 1 for r > ^. At the 
transition point convergence to the scaling form proceeds exceedingly slow, ~ 1/ In t, in 
a fashion reminiscent of critical slowing down in equilibrium phase transitions. 

Several applications come to mind. For example, complex networks could be grown 
according to this model where the nodes are fixed at the outset (corresponding to the A 
agents) and links are connected to the nodes by a proper mix of homogeneous selection 
and preferential attachment. For r > | one could thus create scale- free nets with a 
fixed degree distribution exponent and a fixed number of nodes, and with a tunable 
average connectivity (k) = t/A that grows linearly with time. Wealth distributions 
with a stretched-exponential decay on one side and a power-law decay on the other, 
such as we find for r > |, are regularly observed in various economic settings [HI [9]. 

An intriguing finding concerns the method of rate equations, or mean-field method 
that is often used to obtain the degree distribution of complex networks p, [T0l[m[r2| [T3]. 
Our analysis suggests that this method is only valid when the second-order terms in the 
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Kramers-Moyal expansion of the master equation for the system may be safely neglected. 
Even then the method yields results that scale correctly but that are otherwise only 
qualitatively correct, at least in our case. Perhaps the most important open problem is 
to establish the range of validity of the rate equation approach more rigorously, and to 
find ways to extend it to the cases where it fails. 
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